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Abstract 

The classical Duhamel principle, established nearly 200 years ago by 
Jean-Marie-Constant Duhamel, reduces the Cauchy problem for an in- 
homogeneous partial differential equation to the Cauchy problem for the 
corresponding homogeneous equation. Duhamel's principle is not appli- 
cable in the case of fractional order differential equations. In this paper 
we formulate and prove fractional generalizations of this famous principle 
directly applicable to a wide class of fractional order differential-operator 
equations. 



1 Introduction 

Let AT be a reflexive Banach space and A : V ^ X a, closed linear operator with a 
domain T) C X. In Section 2 we will introduce a Frechet type topological vector 
space Exp^ q(X) (and its dual Exp^, (X*)), where G is an open subset of 
the complex plain C. This space represents a modification of the space of entire 
functions with finite exponential type [G] [H [21] [26] and its abstract versions. 
We also introduce a functional calculus in the form f{A), where / is an analytic 
function defined on G. The function / is called the symbol of the operator f{A). 

The goal of this paper is to generalize Duhamel's principle for the Cauchy 
problem for general inhomogeneous fractional distributed order differential- 
operator equations of the form 

L^[u]= f{a,A)D'^u{t)dA{a) = h{t), t > 0, (1) 







u'-''\0) ^ ipk, fc = 0,...,m-l, (2) 

where ^ € (m — 1, to]; h{t) and (pk, A; = 0, to — 1, are given A-valued vector- 
functions; f{a,A) is a family of operators with the symbol f{a,z) continuous 
in the variable a S [0, and analytic in the variable 2:SG'cC;Aisa finite 
measure defined on [0,/Lt]; and is the operator of fractional differentiation of 
order a in the sense of Caputo-Djrbashian (see, for example, |S1[I1]), i. e. 

D'^qlt) = Jr>^/o ifn-l<a<n, neN, 

\gir^){t)^j^g{t), if« = nG{0}UN. 
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Hereafter the integrals are understood in the sense of Bochner if g{t) is a vector- 
function with values in some topological- vector space for each fixed t. 

The classical Duhamel principle is not applicable in the case of fractional 
order differential equations. Its modification combined with some integral trans- 
formations can reduce the Cauchy problem for inhomogeneous equation to the 
Cauchy problem for homogeneous equation. However, this two step process 
becomes cumbersome for complex equations containing many terms with frac- 
tional operators. In this paper we formulate and prove fractional generalizations 
of Duhamel's principle applicable directly to the Cauchy problem for inhomoge- 
neous fractional order differential-operator equations, which reduce them to the 
Cauchy problem for corresponding homogeneous equations. In the particular 
case of fractional order partial differential equations with a single "fractional" 
term in the equation ([1]), a fractional analog of Duhamel's principle was obtained 

in [Si 137]. 

Fractional order differential equations arc useful and appropriate mathemat- 
ical apparatus for modeling problems with memory, and interest in this subject 
has grown substantially during the last a few decades. For instance, probabil- 
ity density functions of a wide class of non-Gaussian diffusion processes satisfy 
fractional order governing equations with space and time fractional order dif- 
ferential operators (see [U [51 [TSl [Mj and references therein) . Inhomogeneous 
fractional order differential equations appear naturally describing the influence 
of an external force or memory effects. In the study of diffusion processes in 
complex heterogeneous media with several distinct diffusion modes, even with- 
out an external force, the function h{t) embodies memory of the past [22 l [38]. 

There is extensive literature on the Cauchy problem for integer order abstract 
differential-operator equations (see, e.g. [20l [40]). The first order evolution 
equations u (t) ~ Au{t) in the spaces of abstract exponential vector-functions 
of a finite type, Expa{X) (and in more general bornological spaces) were studied 
in [26]. In the case of integer a^, fc = 1, ...,to, the Cauchy problem for pseudo- 
differential and differential-operator equations with analytic symbols or with 
symbols having singularities was studied, for example, in [71[3TJ[52], and multi- 
point value problems in [551 1311 1331 IM]- What concerns fractional order 
differential-operator equations, Kochubei [T7] studied existence and uniqueness 
of a solution to the abstract Cauchy problem D'^uit) = Au{t), u(0) = ito, with 
Caputo-Djrbashian fractional derivative for < a < 1 and a closed operator 
A with a dense domain T>{A) in a Banach space. El-Sayed [11], Bazhlekova [3] 
investigated the Cauchy problem for < a < 2. In the more general case of 
a > 0, Gorenfio et. al. [Tl] studied existence of solutions in Roumieu-Beurling 
and Gevrey classes. Kostin [19] proved correctness of the abstract initial value 
problem (Cauchy type problem) £>"u(i) = Au{t), D"~'^u{0) = ipk, k = I, m, 
for a € (m — l,m), and with the Rieman-Liouville derivative D'^. For more 
information about recent results on the Cauchy problem for abstract fractional 
differential-operator equations, we refer the reader to [31[Tni[IS]; and for a recent 
mathematical treatment of the distributed fractional order differential equations 
to papers [T8 l [23 l [35] . 

The paper is organized as follows. In Section 2 we recall the classic Duhamel 
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principle, and the basic spaces of elements used in this paper. Since we formulate 
a fractional Duhamel principle in the abstract case, we introduce a topological- 
vector space on which the corresponding operators act. In Section 3 wc formu- 
late the main result, namely an abstract fractional analog of Duhamcl's principle 
and discuss some of its applications. 



2 Preliminaries 

2.1 Fractional order derivatives. 

For a function g defined on [0, oo), under some integrability conditions the frac- 
tional integral of order /? with terminal points t and t, is defined as 

rJ^9it) = Y^^j\t-sf"'g{s)ds, 

where r(-) is Euler's gamma-function. Obviously, if /3 = n then ^-J" is the 
n-fold integral of / over the interval [r, i]. By convention, T-J°/(t) = f{t), i.e. 
T- J° coincides with the identity operator. In the notation we do not indicate the 
upper terminal point t, since in the current paper it is always t. 

Further, let m be a positive integer number. We denote by tD", m — 1 < 
a < 1, the fractional derivative of order a in the sense of Riemann-Liouville, 
which is defined as 

g{s)ds 

r(m - a) dt"'^ {t - sY 

and rD+git) = g{t), rD+g{t) = Between this fractional derivative and 

the Caputo-Djrbashian derivative there is the following relationship [12) : 



^DU{t) = — ^ / , ^^ 'Z , m - 1< a < m. 



m — 1 



(fc) ( \ 

rDlgit) = rD:g{t) + ^ - r)^—, t > 0. (3) 

In the particular case of < a < 1 one has 

rD^g{t) ^ rD^9{t)+gir) ^^^^^]'Jy t > 0. (4) 

If .g(r) = 0, then one obtains the equality rD'^+g{t) = T-D"(?(t). Alternative 
representations via the fractional integral are: 

^D^g{t) = —rJ"'-g{t) and ,D^g{t) = ,J'^-^S2, 

We omit the lower terminal point r if r = 0, writing simply Z?" , D" or J". Recall 
that for the Laplace transform of D'^g{t) and D"g{t), where a £ (m — l,m], 
the following formulas are valid P^ : 

fc=0 ^ ' 
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and 



m — 1 



C[D-g]{s)^s'^C[g]{s)^Y.9^'''(^+>' 



.a—l — k 



(6) 



fe=0 



Here C[g]{s) denotes the Laplace transform of g. 
2.2 An operator calculus 

In this section we recall some necessary facts about abstract spaces of analytic 
elements of finite exponential type, and an operator calculus defined on it. See 
for details [32l|33]. 

Let X be a reflexive Banach space with a norm ||w|| , v G X. Let A be a closed 
linear operator with a domain V (A) dense in X and a spectrum (7{A) C C. 
Assume that cf{A) is not empty and is not bounded. 

We will develop an operator calculus f{A) for analytic functions /(A) in an 
open domain G C C If the domain of analyticity of /, G contains cr(^) then 



where u is a, contour in G containing (t{A), and TZ{C,A), ( & C \ (7{A), is the 
resolvent operator of A. However, if / has singular points in the spectrum of A, 
then f{A) can not be defined through the integral ([7]). 

Assume that G is any open set in C not necessarily containing (y{A). Further, 
let < r < +00 and ly < r. Denote by Exp^ ^(X) the set of elements v G 
nk>i'D{A'^) satisfying the inequalities < Ci''^||w|[ for all k = 1,2,..., with 

a constant C > not depending on k. An element v € Exp^ is said to be 

a vector of exponential type v [2^. A sequence of elements w„, n = 1, 2, is 
said to converge to an element vq in Exp^ iff; 

1) All the elements w„ are vectors of exponential type v < r, and 

2) \\vn - foil ^0, n ^ oo. 

Obviously, Exp^ (X) C Exp^ 1/2 (-^)' ^1 < '^2- Let Exp^ ^(X) be the induc- 
tive limit of spaces Exp^ uOQ when v ^ r. For basic notions of topological 
vector spaces including inductive and projective limits we refer the reader to 
P7] . Set Ax ^ A - A/, where A G G, and denote Exp^^_;^(X) = {ua G AT : 
u\ G Exp^^ r(-^)}: with the induced topology. Finally, for arbitrary G C o-{A), 
denote by Exp^ gOQ ^^"^ space whose elements are the locally finite sums of 
elements in Exp^ ^ xiX), A G G, r < dist{X,dG), with the corresponding topol- 
ogy. Namely, any u G Exp^ gO^) ^ representation u = u\ with a finite 
sum. It is clear, that Exp^ gO^) is a subspace of the space of vectors of expo- 
nential type if r < +00, and coincides with it if r = +00. Moreover, Exp^ gO^) 
is an abstract analog of the space \E'g,p(-R^) introduced in [32], where A = — 
G C X= Lp{R^), 1 < p < 00. In the case A = X= L2{R^), the 

corresponding space was studied in [7]. 

Further, let /(A) be an analytic function on G. An arbitrary clement u G 
Exp^ QpC) is represented as a finite sum u = X^agg"^' ^ Exp^^^ ;^(X). 




(7) 



4 



Then for u £ Exp^ gO^) operator f{A) is defined by the formula 

fiA)u - /a(^)"a, (8) 



xeG 

where 



f,(A)u, = Y,^^A-,u,. (9) 

n=a 

In other words, each f\ represents locaUy / in a neighborhood of A G G, and 
for u\ the operator f\{A) is well defined. 

Additionally assume that there exists a one-parameter family of bounded 
invertible operators U\ : X X such that 

AUx - UxA = XUx, A e a{A). (10) 

For example, ii X — L2 = L2{R) and A = — : L2 ^ L2 is the operator of 
differentiation with domain T>{A) = {v ^ L2 Av £ L2}, then for the operator 
Ux '■ v{x) — e*^^f (x) we have 

AUxvix) = -^{e'^'^vix)) = Xe'^'^vix) - te'^^^ = 
ax ax 

XUxv{x) + UxAv{x), 

obtaining (fTOl) . Condition (fTO| indicates a shift of the spectrum of operator A 
to A. This is seen from the relationship A — XI = UxAU^^, which follows from 
([To]) multiplying by U^^ from the right. It follows from the latter that 

(A- A/)" = C/aA"[/-S (11) 

for all n = 1,2, .... 

Let X* denote the conjugate of X, and A* : X* — s- X* be the operator conju- 
gate to A. Further, denote by Exp^, q, (X*) the space of linear continuous func- 
tionals defined on Exp^q(X), with respect to weak convergence. Specifically, 
a sequence G Exp^, q,(X*) converges to an element u* G Exp^. q,(X*) if 
for all V G Exp^ aOQ the convergence < u'^ — u,v >^ holds as to — > cxo. 
For an analytic function /™ defined on G* = {2; G C : z G G}, we define a weak 
extension of f{A) as follows: 

< r{A*)u*,v>^< u\ f{A)v >, V« G Exp^^g(X). 

Lemma 2.1. Let X he a reflexive Banach space and A be a closed operator 
defined on ^{A) C X. Let f be an analytic function defined on a open connected 
set G C C. Then the following mappings are well defined and continuous: 

1. /(A):Exp^^g(X)^Exp^^g(X), 

2. r(A*) : Exp;,. G.(X*) Exp;,. G^X*). 
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Proof. Notice that f{A) maps Exp^(3(X) into itself. Let u £ Expy^Q{X.) 
has a representation u = "^^ux, ua G Exp^^ Then for f{A)u defined in 

([5]), one has the following estimate 

\\AyxiA)u,\\ <f2\I^\\iA-XirAlux\\ < Cz.'=||u,||. (12) 

ri=0 

with some < r. It follows that f\{A)u\ G Exp^^ with the same v, and 
f{A)u G Exp^ gO^)- The estimate (fT2|) also implies continuity of the mapping 
f{A) in the topology of Exp^^ c.(X). 

Now assume that a sequence u* e Exp^. (X*) converges to in the weak 

topology of Exp^. Q, (X*). Then for arbitrary u G Exp^ gO^) '^^ have 

< r{A*)u;^,u >=< KJ{A)u >=< 

where v = f{A)u G Exp^ gO^) '^^'^ the first part of the proof. Hence, 
— > 0, as n — > oo, in the weak topology of Exp^. ■ 

Remark 2.2. It is not hard to see that the above constructions are valid with 
corresponding specifications in the case of operators with discrete spectrum as 
well. Note that in this case the space Exp^ gOQ consists of the root lineals of 
eigenvectors corresponding to the part of a (A) with nonempty intersection with 
G. If the spectrum (j{A) is empty then an additional exploration is required for 
solution spaces to be non-trivial (for details see, \3I)- 

As is shown in [33], the space Exp^ g(-^) invariant with respect to the 
action of an operator f{A) and this operator acts continuously. 

2.3 Two lemmas 

The following two lemmas will be usful in proofs of theorems in Section 3. 

Lemma 2.3. Let h[t) be a continuous differentiable function. Then the equation 
J°'u{t) = h{t), t > 0, where < a < 1, has a unique continuous solution given 
by the formula 

u{t) = Dl-°'h{t), t > 0. (13) 

Lemma [273l is essentially the well-known result on a solution of Abel's integral 
equation of first kind. See [12l [28] for the proof. 

Lemma 2.4. Suppose v(t,T) is a vector-function in a Banach space X, defined 
for all t > T > and k times differentiable with respect to the variable t. Let 
uit) = v{t,T)dT. Then 

-^u(t) = y — [— — — w(i,T)| J + / —rv(t,T)dT. (14) 

dtf" j-'^dP^dt^-^-i Jo df" ' ^ ' 
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Proof. For a fixed t > and small h one can easily verify that 
u(t + h)-u(t) 1 , , , , , /■* ^ , ^ 

^ ' v{t + h,T)dT- v{t,T)dT) 

'o Jo 



' v{t,T)dT + / ^ { —dr 



h 

_ 1 
^ h 

-^r-'C+'-'-Wd.. (15) 

Jt h 

Making use of the mean value theorem (in the integral form), we obtain 



t+h 

v{t,T)dT -V{t,t)\\ < Cl||w(i,T,) --V{t,t)\\, t<T^<t + h, (16) 

'"(' + ")-"<"*-/' ^.,11 <C,W, (IT) 



h Jq dt 

v{t + h)-v{t) 



dT\\ < c^H (18) 

where constants Ci, C2, and C3 do not depend on h. Now, letting /i — s- 0, 
estimates (fT6| -(fT8 | and equation (fT5|) imply the following formula: 

-U{t) = V{t, t) + -V{t, T)dT. (19) 

Formula (|14p follows from (|19p by differentiation repeatedly. ■ 



2.4 Classical Duhamel's principle. 

Duhamel's principle was formulated first for the Cauchy problem for second or 

d_ 

dt 



der linear inhomogeneous differential equations. Let B = B{x, -S^jD^), where 



Dx = ( g§7: be a linear differential operator with coefficients not de- 

pending on t, and containing temporal derivatives of order not higher than 1. 
Consider the Cauchy problem 

-^{t,x)+Bu{t,x)^ f{t,x), t>0,xeR", (20) 

with homogeneous initial conditions 

u(0,x)=0, ^(0,.t)=0. (21) 

Let a sufficiently smooth function v{t, r, x), t > t, t > 0, 2; e i?", be for t > t 
a solution of the homogeneous equation 

— {t,T,x)+Bv{t,T,x)=0, 
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satisfying the following conditions: 



dv 

v{t,T,x)\t=r = 0, —{t,T,x)\t=r = h{T,x). 



Then a solution of the Cauchy problem ((20)) . ((2T|) is given by means of the 
integral 

u{t,x) = / v{t,T,x)dT. (22) 
Jo 

The formulated statement is known as Duhamel's principle, and the integral in 
(|22[) as Duhamel's integral. 

A similar statement is valid in the case of the Cauchy problem with a ho- 
mogeneous initial condition for a first order inhomogencous partial differential 
equation 

du 

— {t,x) + Cu{t,x)^f{t,x), t>0, xGi?", 

where C = C{x,Dx:) is a linear differential operator containing only spatial 
derivatives, and with coefficients not depending on t (see [4]). 



3 Generalizations of Duhamel's principle 

In this section we prove abstract fractional generalizations of Duhamel's princi- 
ple and discuss some of their applications. 



3.1 Duhamel's principle: A = YlT=o^ai, with = k, k = 
1, ...,m. 

Suppose the measure A in ([T|) has the form A = J2T=o : where 6a denotes 
Dirac's delta with mass on a. Suppose also that f{m,A) = /, the identity 
operator. Then the Cauchy problem ([T]), ^ takes the form 

m — 1 

u^"'Ht)+Y,fk{A)u^''Ht) = h{t), t>0, (23) 

uW(0) = fc = 0,...,TO-l. (24) 

The operators fkiA) = f{k,A),k = 0,...,m — 1, are understood in the sense 
of the functional calculus introduced in Section 12.21 In the following theo- 
rem we assume that the vector-functions U{t,T) and h(t) are Exp^^ q(X)-, or 

Exp^. (X*)-valued. In this abstract case Duhamel's principle is formulated 
as follows. 
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Theorem 3.1. Let a vector-function U{t, r) for all t : < t < t be a solution 
of the Cauchy problem for a homogeneous equation 



f)mTj ™ 1 rl^TJ 

{t^r)+Y,f,{A)-^{t,T)=Q, t>r, 



k=0 



(t,T)\t=T+a = 0, fc = 0, ...,m-2, 



Qra-ljj 

{t,T)\t=r+0=h{T), 



(25) 

(26) 
(27) 



where h{t) is a continuous vector-function. Then a solution of the Cauchy prob- 
lem for the inhomogeneous equation 

rn— 1 



(t) + ^/,(A)uW(t) = Mt), 



fe=0 



^^='(0) =0, k^O,...,m-l. 
is represented via Duhamel's integral 

u{t) = [ U{t,T)dT. 



(28) 
(29) 

(30) 



Proof. Obviously u(0) — 0. Further, for the first order derivative of w(i), 
using Lemma \2M one has 



^ii) = m,t)+l ^it,r)dr, 



By virtue of (l26l) the latter implies that ^(0) ~ 0. Further, differentiating. 



( 

'dt^^*'~ df^ 



— {t,t) + 



dt^ 
(0) = 0, fc = 2,...,m- 1 



it,T)dT, 



which due to condition (p6|) implies that 

d'^u , 
IF 

Therefore, the function u{t) in (pO)) satisfies initial conditions ((29)) . Moreover, 
substituting ([50)1 to ([^5]). and taking into account (P7)) . we have 



^(™)(i)+^/,(A)uW(t)^ / U{t,T)dT+J2f,{A)— U{t,T)d7 



■m — 1 



A:=0 

gm-iu 
/i(i) 



fc=0 



-(i,r)dr+^/fc(A) 



at" 



5r 

(t,r)+^/,(A)^(t,r; 



fc=o 



k=0 



df^ 
dr = /i(t) 



{t,T)dT 



Hence, u{t) in (|30l) satisfies equation (|28l) as well. 
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Remark 3.2. It is not hard to see that Theorem \3. 1\ holds with generic closed 
operators Bk (with dense domain 'D{Bk) and commuting with -^j instead of 

fk{A). In this case we assume that hit) £ X and ^ gj'f''^'' S V{Bk), k = 0, m— 



3.2 Fractional Duhamel's principle: A = 6^ + X with fi G 

(m — 1, m]. 

Let A — 6fj, + X, where /.i is a number such that m — I < fi < m, and A is a finite 
measure with suppX C [0,m — 1]. Consider the operator 

prn—l 

rL^'''^\u]{t)= rD>:u{t)+ f{a,A)rDMt)Hda), (31) 



acting on m-times difFerentiable vector-functions u{t), t > r > 0. If r = 0, then 
instead of oi^^'^'N(i) we write 

Theorem 3.3. Suppose that V{t,T), t > t > 0, is a solution of the Cauchy 
problem for the homogeneous equation 

.L('^'^)[y(.,r)](i) = 0, t>T, (32) 

-^{t,T)\t=r+0^0, fc = 0,...,TO-2, (33) 

^-^it,r)U=r+o=Dr''h{r), (34) 
where h{t) is a given vector-function. Then Duhamel's integral 

U{t) = [ V{t,T)dT (35) 

Jo 

solves the Cauchy problem for the inhomogeneous equation 

L(^-^)[u](t) = h{t), t > 0, (36) 

with the homogeneous Cauchy conditions 

u('^')(0)=0, fc = 0,...,m-l. (37) 

Proof. First notice that since to — 1 < < m, and therefore < m — /i < 1, 
due to Lemma [2.31 the equation J™'^'^g{t) = h{t) has a unique solution 

g{t) = D'r'Kt). (38) 

Let V{t, t) as a function of the variable t be a solution to Cauchy problem 
(|32|) - (p4| for any fixed r. We verify that u{t) — V{t,T)dT satisfies equation 
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and conditions ([57)) . Splitting the interval (0,m — 1] into subintervals 
(0, 1], (m — 2, m — 1], we have 

x(^'^)M(t) = + V / f{a,A)Dy{t)X{da). (39) 

fc=i ^fe-i 

If a G (fc — 1, fc], where fc = 1, m — 1 using the definition of D", we have 



Lemma [2.41 and conditions ([33]) imply that 
''" lo Jo 



Jk l-s rs Qk 

j V{s, r)dr = J ^V{s, r)dr, fc = l,...,m-l. (41) 



Hence, for aG(fc — l,fc],fc=:l,...,TO — 1, 

D^uit) = ^ j\t - sr-^^Vis, r)dsdr. (42) 
Again due to Lemma and condition 



Therefore the first term on the right hand side of ([39)) takes the form 



r(m 



1 (7™ 



r(TO - ^l) Jo 

Further, put g{t) ^ Dl'^^h{t). Then by virtue of ([38)). 



(t-.)™-''-l(i?r^M.) + I —V{s,r)dT)ds. 

(44) 



1 



(i „ s)™-^-iZ?"^^/i(s) = J"-''.g(i) = /i(t). (45) 



r(m - ^() Jo 

Now equations ([51)), ([H]), OH)), and 015)) imply that 



L(^^^'M(0 = Mt) + -F7Z-^ / (^-^r-^"' / ■T—V{s,r)drds 



r(m - /i) Jo Jo ds 

m-l »fe »t 



(46) 



k=l 
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Changing the order of integration (Fubini is aUowed) in we get 



t l-t 1 an 

1^ „\m-/i-l <^ 



™-l ("fc ft 1 ofe 

= h{t)+ rD^^-V{t,T)dT + / ,7^?l/(t,T)A(da)dT 
Jo Jo Jo 



Finally, using the relations (|4T|) it is not hard to verify that u{t) in (|35|) satisfies 
initial conditions ([57)) as well. ■ 

If the vector-function /i satisfies the additional condition /i(0) = then 
condition ([M)) in accordance with the relationship can be replaced by 

Qrn-ly 

——it^r)\,=^^Df~^h{r), 

with the Caputo-Djrbashian derivative D"^^^ of order m — ji. As a consequence 
the formulation of the fractional Duhamel's principle takes the form: 

Theorem 3.4. Suppose that for all t : < t < t a function V{t,T), is a 
solution to the Cauchy problem for the homogeneous equation 

rL(^'^'^[V{;Tm=0, t>T, 

-^{t,T)\t=r+0 ^ 0, fc = 0,...,TO-2, 
gm-ly 

^—^it,r)U^.^o = DT-^'hir), 
where h{t) is a given vector-function such that h{0) = 0. Then 

u{t) = / V{t,T)dT 

Jo 

is a solution of the Cauchy problem for the inhomogeneous equation 

L^i'^^\u]{t)^h[t), i>0, 
with the homogeneous Cauchy conditions 

u('=)(0) = 0, fc = 0,...,m-l. 



Remark 3.5. 1. Lemma \2.S\ can be extended to absolutely continuous func- 
tions h{t) with an appropriate meaning of solution in equation (|13[) (see, 
\28f ) . It is also known \28f that the fractional derivative D^^h{t), k — 
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1 < /i < fc, fc = l,...,m, exists a.e., if h(t) is an absolutely continuous 
function on [0; T] for any T > 0. These two facts imply that general- 
ized Duhamel's principles proved above hold true for absolutely continuous 
functions h{t). 



2. In Theorems \3.3\ and \3.4\ we assumed that f[^,A) is the identity operator 
(see equation (pij) ). In the general case, with appropriate selection of G 
we can assume that the inverse operator [f[ii,,A)\~^ exists. Then with the 
condition 



Qtm-l \t=T 



= [/(M,A)]-ii?r^Mr) 
instead of (|34p . Theorems \3.3\ and \3.4\ remain valid. 

3.3 Fractional Duhamels's principle with Riemann-Liuoville 
derivative 

The operator t-L^ in Theorem 13. 31 is defined via the fractional derivative in the 
sense of Caputo-Djrbashian. A fractional generalization of Duhamel's principle 
is also possible when this operator is defined via the Riemann-Liouville fractional 
derivative. In this section we briefly discuss this important case proving the 
corresponding theorem in the simple particular case TL[u]{t) = rD+u(t) + 
Bu{t), where < a < 1, and i? is a closed operator with a domain V{B) dense 
in X. The general case can be treated in a similar manner. 

Theorem 3.6. Suppose that V{t,T), t > t > 0, is a solution of the Cauchy 
type problem for the homogeneous equation 

rD^V{t,T) + BV{t,T) =0, t>T, (47) 
rJ'-''V{t,T)\^^^^^h{T), (48) 

where h{T), r > 0, is a continuous vector-function. Then Duhamel's integral 

U{t) = / V{t,T)dT (49) 

Jo 

solves the Cauchy type problem for the inhomogeneous equation 

D%u{t) + Bu{t) = h{t), t > 0, (50) 

with the homogeneous initial condition J^^"u(0) = 0. 

Proof. Let V{t, r) satisfy the conditions of the theorem. Then for Duhamel's 
integral , by virtue of Lemma 12.41 we have 

1 d /■* CVis,T)ds /■* 

= ^ I rJ^-"V{t,T)dT+ [ BV{t,T)dT 

ax Jo Jo 

= ,j'-"V{t,T)\^_^ + f [rDlV{t,T)+BVit,T)]dT = h{t). (51) 

Jo 
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On the other hand, changing the order of integration and using the mean value 
theorem, we obtain 

= II r rJ'-''V{t,T)dT\\<t\\rj'-''V{t,T)\\. (52) 

Jo 

Condition (^5)) implies that limt^o+ TJ^^°'V{t,T) = h{0) in the norm of X. It 
follows from ([5^ that limt_^o+ J^^°'u{t) = in the norm of X. ■ 

3.4 Applications. Existence and uniqueness theorems 

Theorems 3.1 and 3.3 lead to generalization of the existence and uniqueness 
results obtained in papers [131 [30] for the abstract Cauchy problems. Let 
be the distributed fractional order abstract differential-operator defined in (|3ip 
with T = 0, and via the characteristic function 

/•m— 1 

A{s,z) = s^+ f{a,z)s"dX, 
Jo 

where fi G (m—1, to], A is a finite measure with supp Ac [0, m — 1], and f{a, z) is 
a fimction continuous in a and analytic in z G G C C Denote by v{s) ~ C[v]{s) 
the Laplace transform of a vector- function w(t), namely 

POO 

C[v]{s)^ / e^-'^vitjdt, s>so, 
Jo 

where sq is a real number. It is not hard to verify that if v{t) G ExpA^ai^) 
for each t > and satisfies the condition ||w(i)|| < Ce^^^t > 0, with some 
constants C > and 7, then v{s) exists and 

s — 7 

implying £1(5) G Expa,g{X) for each fixed s > 7. The lemma below gives a 
formal representation formula for a solution of the general abstract Cauchy 
problem 

L^[u]{t) ^ h{t), t > 0, (53) 
it('=)(0+) = ipk, fc = 0,...,TO-l. (54) 

Lemma 3.7. Let cp{t,z) = ^^^[-^^jj^Kt), z lE G G C, where Ct^ stands for 
the inverse Laplace transform, and 

pm — l 

Sk{t,z) = Cf,^k-i{t,z) + J f{a,z)ca-k-i{t,z)X{da). 

Then Uk(t) — Skit, A)ipk solves the Cauchy problem 

L[u] = 0, u\Q) ^ d.j^k'P], j = 0,...,m-l, 
where Sj,k = 1 if j ~ k, and 5j^k ~ 0, if j ^ k. 
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Corollary 3.8. Let Sk{t, A), fc = 0, m — 1, be the collection of solution oper- 
ators with the symbols Sk{t,z) defined in Lemma \3. 7[ Then the solution of the 
Cauchy problem 

L[u]=0, u^(0) =v3j , j = 0, ...,m-l. (55) 
is given by the following representation formula 

m — 1 

u{t)=Y,Sk{t,A)^k- (56) 

A:=0 

Proof of Lamma. Applying formula ^ we have 

m— 1 



C[L[u]]{s) = sf'uis) ~ J2 u\0)s''-'-^ + 

i=Q 

/ /(a,A)(s"u(s)-^u-'"(0)s"-^-i)A(da) =0. 

fc = l "''=-1 7=0 



Due to the initial conditions m-' (O) = Sj^k'Pj, j — 0, m — 1, the latter reduces 
to 



' m — 1 

A(s, z)uis) = ( s^-'-^ + J /(a, z)s"-'=-iA(da) 
Now it is easy to see that the solution in this case is represented Uk = Sk{t, A)ipk. 



Remark 3.9. 1. Corollary \3.8\ can easily be extended to the operator ^-L in 
(|55p as well with the initial conditions u\t) = ipj, maintaining the shift 
t = t~T; 

2. A particular case of Lemma \37)\ when A = J2T=o , k — 1 < au < k, is 
proved in U3^ . 

Further, denote by C('")[t > 0; Expam{X)] and by AC[t > 0; ExpA,GiX)] 
the space of m— times continuously difFerentiable functions and the space of 
absolutely continuous functions on (0; +oo) with values ranging in the space 
Expa,g{X), respectively A vector-function u{t) e C(™^[t > 0; ExpA^aiX)] 
n C(™-i)[t > 0; Expa,g{X)] is called a solution of the problem 1^, if 
it satisfies the equation ([53]) and the initial conditions (|54| in the topology of 
ExpA^ciX). 

Theorem 13.31 and Corollarv 13.81 implv the following results. 

Theorem 3.10. Let ipk £ Expa,g{X), fc = 0, ... , to - 1, h{t) e AC[0 < t < 
T; Expa,g{X)] for any T > 0, and D'^-'^-^hit) eC[0<t<T; Expa,g{X)]. 
Then the Cauchy problem i53\) . ^54^ has a unique solution. This solution is 
given by 

m-l -t 

n{t) = Sk{t,A)^k + / - T,A)D'^'^'h{T)dT. (57) 

fe=o -^0 
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Proof. Wc split the Cauchy problem ([5^ . ([5^ into two Cauchy problems 

L^[u]{t) 0, t > 0, (58) 

u(^-)(0+) = ^fc, fc = 0,...,77i-l, (59) 

and 

L^[v]{t)^h{t),t>0, (60) 

vW(0+) 0, fc = 0, ...,m- 1. (61) 

Due to Corollarv l3.8l tlie unique solution to ([55)) . ([5^ is given by 

m— 1 

u{t)=J2Sk{t,A)^k- (62) 

fe=0 

For the Cauchy problem (|60p . (|6T|) . in accordance with the fractional Duhamel's 
principle (Theorem 13. 3p . it suffices to solve the Cauchy problem for the homo- 
geneous equation: 

rL^[V{t,T)]{t) ^0,t>T, (63) 



0, fc = 0,...,m-2, (64) 



|t=r + 



^ ' ' D'^-^hir). (65) 



at" 



t = T + 



The solution of this problem, again using Corollarv l3.8l (with the note in Remark 
13. 9|) . has the representation 

V{t,T) = 5™_i(t-T,^)Z?™-''/i(r). (66) 

Thus, Duhamel's integral of the latter and representation lead to formula 
(|57|) . The uniqueness of a solution also follows from the obtained representation 
(EZD (see [33). ■ 

The duality immediately implies the following theorem. 

Theorem 3.11. Let ^pl e Exp'^, (^,{X*), k = 0, ... , m - 1, h*{t) e AC[t < 
t < T- Exp\, (.,{X*)] and Dl'~""^h*(t) e C[t < t < T; Exp'^, c!,{X*)]. 
Assume also that Expa,g{X) dense in X. Then the Cauchy problem \53\] . 
{54-^ (with A switched to A*) is meaningful and has a unique weak solution. 
This solution is given by 

U*{t)=Y,Sk{t,A*M+ / S,n-l{t-T,A*)D^-''h*{T)dT. 

fc=0 -^0 
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Assume that Expa,g{X) is densely embedded into X. Besides, let the solu- 
tion operators Sk{t, A) for each k = 0, m — 1, satisfy the estimates 

<C||^||, V<e[0,T], (67) 

where e Expa,g{X), and C > does not depend on ^p. Then there exists a 
unique closure Sk{t) to X of the operator S'fc(<, A) which satisfies the estimate 
||<5A;(t)M|j < IImII for all u E X. Using the standard technique of closure (see 
we can prove the following theorem. 

Theorem 3.12. Let tpk E X. k ^ 0, ... , m - 1, h{t) G AC[0 < t < T; X] 
for any T > 0, and D'^^°''^h{t) e C[0 < i < T; X]. Further let Expa,g{X) 
be densely embedded into X, and the estimates \61^ hold for solution operators 
Sk(t,A),k ~ 0, ...,m — 1. Then the Cauchy problem i53\} . |5^p has a unique 
solution u{t) G C™[0 <t<T;X]. This solution is given by 

m— 1 
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